Rectification and nonlinear transport in chaotic dots and rings 
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We investigate the nonlinear current-voltage characteristic of mesoscopic conductors and the cur- 
rent generated through rectification of an alternating external bias. To leading order in applied 
voltages both the nonlinear and the rectified current are quadratic. This current response can be 
described in terms of second order conductance coefficients and for a generic mesoscopic conductor 
they fluctuate randomly from sample to sample. Due to Coulomb interactions the symmetry of 
transport under magnetic field inversion is broken in a two-terminal setup. Therefore, we consider 
both the symmetric and antisymmetric nonlinear conductances separately. We treat interactions 
self-consistently taking into account nearby gates. 

The nonlinear current is determined by different combinations of second order conductances de- 
pending on the way external voltages are varied away from an equilibrium reference point (bias 
mode). We discuss the role of the bias mode and circuit asymmetry in recent experiments. In 
a photovoltaic experiment the alternating perturbations are rectified, and the fluctuations of the 
nonlinear conductance are shown to decrease with frequency. Their asymptotical behavior strongly 
depends on the bias mode and in general the antisymmetric conductance is suppressed stronger then 
the symmetric conductance. 

We next investigate nonlinear transport and rectification in chaotic rings. To this extent we 
develop a model which combines a chaotic quantum dot and a ballistic arm to enclose an Aharonov- 
Bohm flux. In the linear two-probe conductance the phase of the Aharonov-Bohm oscillation is 
pinned while in nonlinear transport phase rigidity is lost. We discuss the shape of the mesoscopic 
distribution of the phase and determine the phase fluctuations. 

PACS numbers: 73.23.-b, 73.21.La,73.40.Ei, 73.50.Fq 
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I. INTRODUCTION 

A large part of modern physics is devoted to nonlin- 
ear classical and quantum phenomena in various systems. 
Such effects as the generation of the second harmonic 
or optical rectification are known from classical physics, 
while quantum electron pumping through a small sample 
due to interference of wave functions is a quantum nonlin- 
ear effect. Experiments on nonlinear electrical transport 
often combine classical and quantum contributions. A 
macroscopic sample without inversion center 1 exhibits a 
current- voltage characteristic which with increasing volt- 
age departs from linearity due to terms proportional to 
the square of the applied voltage. If now an oscillating 
(AC) voltage is applied, a zero-frequency current (DC) is 
generated. 

If the sample is sufficiently small, quantum effects can 
appear due to the wave nature of electrons. The uncon- 
trollable distribution of impurities or small variations in 
the shape of the sample result in quantum contributions 
to the DC which are random. For a mesoscopic conduc- 
tor with terminals a, (3, ... we can describe the quadratic 
current response in terms of second order conductances 
Gaf3-y They relate voltages Vp. u applied at contacts or 
neighboring gates (3 at frequency w to the current at zero 
frequency at contact a, 
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(1) 



The second order conductances include in detail the role 
of the shape and the nearby conductors (gates). They 



depend on external parameters like the frequency of the 
perturbation, temperature, magnetic field or the connec- 
tion of the sample to the environment. 

We concentrate here on the quantum properties of 
nonlinear conductance through coherent chaotic samples. 
Chaos could result from the presence of impurities (dis- 
order) or random scattering at the boundaries (ballistic 
billiard). Due to electronic interference the sign of this ef- 
fect is generically random even for samples of macroscop- 
ically similar shape. 2-4 When averaged over an ensemble, 
the second order conductances vanish. As a consequence, 
for a fully chaotic sample there is no classical contribu- 
tion to the DC and the nonlinear response is the result 
of the sample-specific quantum fluctuations. 

Interestingly enough, from a fundamental point of view 
these fluctuations of nonlinear conductance are sensitive 
to the presence of Coulomb interactions and magnetic 
field. While interactions strongly affect the fluctuations' 
amplitude, their sign is easily changed by a small vari- 
ation of magnetic flux <E>, similarly to universal conduc- 
tance fluctuations (UCF) in linear transport. More im- 
portantly, without interactions the current (1) through a 
two-terminal sample is a symmetric function of magnetic 
field, just like linear conductance. However, the idea 
that Coulomb interactions are responsible for magnetic- 
field asymmetry in nonlinear current was recently pro- 
posed theoretically 5,6 and demonstrated experimentally 
in different mesoscopic systems. 7-12 (Various aspects of 
nonlinear quantum 13-16 and classical 17 charge and spin 
transport 18 have been discussed later on.) It is useful 
to consider (anti) symmetric second order conductance 



Ga, Qs defined as 

(Gs(t>)\ h d* / /(<&) ± /(-<&) \ 
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where V is a combination of voltages at the gates and 
contacts varied in the experiment and v s accounts for the 
spin degeneracy. We emphasize that, depending on the 
way voltages are varied, experiments probe different lin- 
ear combinations of second order conductance elements 
Gafj-y of Eq. (1). From now on we will simply call Gs, Ga 
conductances and if no confusion is possible leave out the 
expression " second order" . 

In the presence of a DC perturbation the meso- 
scopic averages of antisymmetric 5 ' 6 and symmetric 14 ' 16 
conductances vanish, and it is their sample-to-sample 
fluctuations that arc measured. Experiments are usu- 
ally performed for strongly interacting samples and the 
magnetic- field components Gs , Ga allow one to evaluate 
the strength of interactions. 10 ' 11 In previous theoretical 
works on nonlinear transport through chaotic dots several 
important issues have been discussed using Random Ma- 
trix Theory (RMT). 5 ' 14 ' 16 Sanchez and Biittiker 5 found 
the fluctuations of Q a in a dot with arbitrary interaction 
strength at zero temperature and broken time-reversal 
symmetry due to magnetic field. Polianski and Biittiker 
considered the statistics of both Q a and Gs for arbitrary 
flux $, the temperature T, and the dephasing rate. 14 
The fluctuations of relative asymmetry A = G a /G s and 
the role of the contact asymmetry on this quantity were 
discussed in Ref. 16. The results of RMT approach were 
compared with experimental data of Zumbiihl et al. 10 
and Angers et al. 11 

Previously we considered statistics of Ga,G s for the 
dots where only one DC voltage was varied. However, 
to avoid parasitic circuit effects some experiments are 
performed varying several voltages simultaneously. Sur- 
prisingly, the importance of the chosen combination of 
varied voltages (bias mode) was not addressed before in 
the literature. It turns out that an experiment where 
only one of the voltages is varied 8,11,19 or two voltages 
are asymmetrically shifted 9,10 measure different combi- 
nations of nonlinear conductances Ga/3-y- For example, 
in a weakly interacting dot in the first mode we found 
that G s 3> Ga, 16 but in the second bias mode the fluctu- 
ations of nonlinear current are strongly reduced, so that 

Gs ~ Ga- 
it is also important to generalize the previous treat- 
ment of the nonlinear current to mesoscopic systems bi- 
ased by an AC-voltage at finite frequency. The result- 
ing DC is sometimes called "photovoltaic current". We 
expect that in such mesoscopic AC/DC converters the 
interactions lead to significant magnetic field-asymmetry 
in the DC-signal. The rectification effect of mesoscopic 
diffusive metallic microjunctions was theoretically con- 
sidered by Falko and Khmelnitskii 20 assuming that elec- 
trons do not interact. Therefore, a magnetic-field asym- 
metry was not predicted and was also not observed in 
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FIG. 1: Top: Quantum pumping sources include oscillat- 
ing voltage Vp(w) at the locally applied gate, which slightly 
changes the shape of the dot (shown dashed), or microwave 
antenna emitting photons with energy huj into the dot. Bot- 
tom: Rectification sources include external bias Vi,2(w), top 
gate voltage Vb(oi) with capacitance C g o, and parasitic cou- 
pling of V p ((j) due to stray capacitances C s tray- Microwave 
antenna can emit photons to the contacts and lead not only 
to photon-assisted AC transport but also to a rectified DC. 



subsequent experiments. 21-25 The fact that the interac- 
tions induce a magnetic field-asymmetry of the photo- 
voltaic current when the size of the sample is strongly 
reduced was recently demonstrated in Aharonov-Bohm 
rings by Angers et al. 12 

However, it turns out that for an AC perturba- 
tion another quantum interference phenomenon, also 
quadratic in voltage, random in sign and magnetic field- 
asymmetric, contributes to the DC. Due to internal AC- 
perturbations of the sample, the energy levels are ran- 
domly shifted and a phenomenon commonly referred to 
as "quantum pumping" 26 ' 27 appears. Brouwer demon- 
strated that two voltages applied out of phase generate 
pumped current linear in frequency, while a single volt- 
age pumps current quadratic in frequency w. 26 Although 
theory usually considers small (adiabatic) frequencies, a 
photovoltaic current could be induced by voltages applied 
at arbitrary frequency. At small uj the pumping contri- 
bution vanishes and only the rectification effect survives. 
In contrast, it is not clear what the ratio of pumping cur- 
rent to rectification current is at large w. To distinguish 
between different mechanisms it is therefore important 
to consider rectification in a wide range of frequencies in 
detail. 

We point here to a crucial difference between recti- 
fication and pumping contributions to the photovoltaic 
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effect. Rectification results from external perturbations 
or the perturbations that can be reduced to the exte- 
rior by a gauge transformation. Typical examples are 
external AC-bias, or gate voltage which shifts all levels 
uniformly, 28 or a bias induced by parasitic (stray) capac- 
itance which connects sources of possible internal pertur- 
bations to macroscopic reservoirs, 29 see the bottom panel 
in Fig. 1 . Pumping, on the other hand, is due to internal 
perturbations like those of a microwave antenna 30 or a 
locally applied gate voltage, 26 see the top panel in Fig. 
1. Internal and external sources affect the Schrodingcr 
equation and its boundary conditions, respectively. In 
experiment pumping and rectification, often considered 
together under the name of photovoltaic effect, 21 ~ 25 - 31 are 
hard to distinguish. 

Can one clearly separate quantum pumping from rec- 
tification effects? To distinguish them it was proposed 
to use magnetic field asymmetry of DC as a signature 
of a true quantum pump effect. In Refs. 29 and 32 rec- 
tification by (non-interacting) quantum dot was due to 
stray capacitances of reservoirs with pumping sources. 
The rectified current was found to be symmetric with 
respect to $ — ► — <E>. 29 While such field-symmetric recti- 
fication dominated in the experiments of Switkes et al. 33 
and DiCarlo et al. 32 at MHz frequencies, an asymme- 
try $ — ■> — $ observed at larger GHz frequencies seemed 
to signify a quantum pump effect. 32 It was noted that 
the Coulomb interactions treated self-consistently do not 
lead to any drastic changes in the mesoscopic distribu- 
tion of a pumped current. 26 Probably, that is why the 
effect of interactions on the rectification have not been 
considered yet, even though the Coulomb interaction in 
such dots is known to be strong. 10 

However, as it turned out later, Coulomb interactions 
are responsible for magnetic-field asymmetry in nonlinear 
transport through quantum dots. 5 Similarly this could 
be expected for rectification as well. Then the mag- 
netic field asymmetry alone can not safely distinguish 
pumping from rectification. Therefore we thoroughly ex- 
amine the frequency dependence of the magnetic-field 
(anti)symmetric conductances Q a ,G s - Here we neglect 
any quantum pumping effects and their interference with 
rectification. 34,35 While the role of Coulomb interactions 
and the full frequency dependence in quantum pumping 
are yet to be explored, here we answer two important 
questions concerning a competing mechanism, rectifica- 
tion: (1) In the DC limit u> — > for a strongly interacting 
quantum dot Q a and Q s arc of the same order. Is this 
also the case at finite frequencies? (2) How are the ex- 
perimental data affected by the bias mode for alternating 
voltages? 

A number of very recent experiments on nonlinear DC 
transport 9 ' 11 and AC rectification 12 have used submicron 
ring-shaped samples with a relatively large aspect ratio. 
In this work we develop a model of a ring which includes 
chaotic dynamics due to possible roughness of its bound- 
ary and/or the presence of impurities. Similarly to quan- 
tum dots, the two-terminal nonlinear conductance of such 



a ring is field-asymmetric because field-asymmetry exists 
in each arm. In particular, this leads to deviations of 
the phase in AB oscillations from O(mod) it which char- 
acterizes linear conductance obeying Onsager symmetry 
relations. Experiments find that the amplitude and phase 
of AB oscillations exhibit rather curious properties. For 
example, the DC experiment of Leturcq et al. 9 finds 
that during many AB oscillations with period hc/e the 
phase is well-defined. The experiment demonstrates that 
a nearby gate can vary the phase of the AB oscillations 
over the full circle. The amplitude of the second har- 
monic hc/2e is strongly suppressed. On the other hand, 
the DC experiment 11 and AC experiment 12 of Angers 
et al. find that the phase can be defined only for few os- 
cillations at low magnetic fields. For high frequencies, the 
phase fluctuates strongly as function of frequency. Both 
in the nonlinear and the rectified current the amplitude 
of the second harmonic hc/2e in AB oscillations is al- 
ways comparable with the first harmonic hc/e. This is 
in contrast with the experiments in Ref. 9. Although we 
do not fully address all these questions here, our model 
of a chaotic ring allows us to consider them at least on a 
qualitative level. 



II. PRINCIPAL RESULTS 

To introduce the reader to the problem of nonlinear 
transport in Sec. IV we first qualitatively discuss the 
Coulomb interaction effect in the simplest DC problem. 
In reality the statistical properties of conductances G a p-y 
in Eq. (1) are sensitive to electronic interference but 
to assess the role of Coulomb interactions we can con- 
sider a specific sample. In contrast to linear transport, it 
turns out that the nonlinear current strongly depends on 
the way voltages at the contacts and/or nearby conduc- 
tors are varied from their equilibrium values (bias mode). 
For example, we find that the experiments when only one 
voltage at the contact is varied 8,11 ' 19 or when two contact 
voltages are shifted oppositely 9 ' 10 measure different non- 
linear currents. Indeed, for a current I({Vi}), bilinear in 
voltages, its second derivative should depend on the cho- 
sen direction in the space of voltages {Vi}. Interestingly, 
a sample with weak interactions is very sensitive to the 
choice of the bias mode, which we attribute to the strong 
effect of capacitive coupling of the sample with nearby 
conductors. 

To make our arguments quantitative and consider the 
role of magnetic flux $ for a quantum dot which is (gen- 
erally) AC-biased at arbitrary frequency w, in Sec. V we 
take electronic interference into account. Having done 
that, we illustrate the interplay between interactions and 
interference on several important examples. First, we 
consider nonlinear transport due to a constant applied 
voltage and then consider rectification of AC voltages. 

For a two-terminal dot, in a generally asymmetric cir- 
cuit (capacitive couplings included), in Sec. VA we find 
the statistics of (anti) symmetric conductances Q ai Q s de- 
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fined in Eq. (2). Both Q a and Q s vanish on average. 
Quantum fluctuations of Q g strongly depend on the inter- 
action strength, circuit asymmetry and bias mode. This 
is in accordance with our qualitative picture. On the 
other hand, the antisymmetric component Q a depends 
only on interactions. Our arguments agree with recent 
experiments in quantum dots: 10 ' 19 depending on the bias 
mode different features of the nonlinear conductance ten- 
sor are probed. The fluctuations of nonlinear current can 
be minimized or maximized (on average), which becomes 
important for weakly interacting electrons. Curiously, 
for symmetric coupling (transmission and capacitance) 
of contacts and dot the bias mode in which the voltages 
at the contacts are changed in opposite directions gener- 
ally minimizes fluctuations of Q s . Consequently, such a 
mode is more advantageous for the observation of Q a or 
a cleaner linear signal. Near the end of Sec. V A we also 
demonstrate how to take into account possible classical 
circuit-induced asymmetry 19 due to the finite classical 
resistance of the wires. 

In Sec. VB we present results elucidating the role 
of interaction in rectification through two-terminal dots. 
Usually there are two important time-scales: the dwell 
time Td an electron spends inside the dot and the charge 
relaxation time trc < T d of the dot. For a given geom- 
etry, the dwell time depends on the coupling of the dot 
with reservoirs, but the charge relaxation time is also 
sensitive to the interaction strength. We have trc <C t<j 
for strong interactions and trc = T d in the weak inter- 
action limit. Our results for fluctuations of Q a (u>), G s (ll>) 
are obtained for arbitrary frequency u>. Although the 
fluctuations of both Q a {u) and Q s (uj) monotonically de- 
crease when uo — > oo, as functions of frequency u they 
behave differently. At nonadiabatic frequencies uot<\ >• 1 
the nearby gate short-circuits currents. This effect is even 
in magnetic field and thus affects only Q s . As a result, 
for a high-frequency voltage the asymptotes of Q a and 
Q s are generally different and strongly depend on the 
bias mode. Since the regime of parameters is quite re- 
alistic, we expect that the predicted difference of Q s (oj) 
and G a {u) should be experimentally observable. In the 
noninteracting limit our results qualitatively agree with 
those in diffusive metallic junctions. 

Our model of a ring consisting of a chaotic dot with 
a ballistic arm which encloses an AB flux is presented 
in Section VI. Although it is impossible to find the full 
mesoscopic distribution of the AB phase 5, its shape can 
be discussed qualitatively. Since tan 5 is similar to the 
asymmetry parameter A — Ga/Gs in quantum dots, its 
distribution can become very wide for a particular choice 
of the bias mode. On average (<5(mod)7r) = in our 
model, and we find the dependence of the fluctuations 
of S on temperature, interactions, and number of chan- 
nels of the contacts and the arm. Our treatment allows 
a straightforward generalization to treat AC voltages ap- 
plied to the ring. The technical calculations are presented 
in the Appendix. 




FIG. 2: (Left) Rectified current is measured through a co- 
herent quantum dot biased by voltages with (AC) amplitude 
Vi tlo ,i = 1,2 at reservoirs connected by iV, ballistic channels 
and capacitances d and by voltages V 9 i^ applied at addi- 
tional gates with capacitances C 3 i. Transport through the 
dot is sensitive to the total magnetic flux $ through the area 
of the dot. (Center and right) Forward and reverse connec- 
tion of Ref. 19 exchange voltages at the contacts and classical 
resistors 7-1,2. 



III. MODEL 

The 2D quantum dot, see the left panel in Fig. 2, is 
biased with several voltages {Vi} at M ballistic quan- 
tum point contacts (QPCs) with Ni,i = 1, ...,M orbital 
channels. The reservoirs can be capacitively coupled to 
the dot via capacitances Cj. An additional set of volt- 
ages {V g i} is applied to (several) gates with capacitances 
C g i. All perturbations are assumed to be at the same 
frequency ui, which is not necessarily small (adiabatic). 

The dot is in the universal regime, 36 when the Thou- 
less energy i?Th = ^/^crg is large. The dots with area 
A = ttL 2 (taken circular) are either diffusive with mean 
free path ! < L, or ballistic, with I 3> L and chaotic 
classical dynamics (in the latter case the substitution 
I —> ttL/4 should be used). The mean level spacing 
(per spin direction) A = 2irh 2 / (m*A) and the total num- 
ber of ballistic channels N together define the dwell time 
r d = h/(NA) > r crg . We also require that eV < iVA 
when we can treat the nonlinearity only to (eV) 2 . Scat- 
tering is spin-independent and this spin degeneracy is 
accounted for by the coefficient v s . 

The noninteracting electrons are treated using the 
scattering matrix approach and Random Matrix Theory 
(RMT) for the energy-dependent scattering matrix 5(e). 
For details we refer the reader to reviews. 36 ' 37 In this 
approach the fundamental property of a dot is its scat- 
tering matrix S distributed over circular ensembles of 
proper symmetry, see Ref. 36 (An alternative method is 
the Hamiltonian approach based on the properties of the 
dot's Hamiltonian H taken from a Gaussian Ensemble. 37 ) 
Transport properties of chaotic dots in RMT for matri- 
ces S or TL are usually expressed in terms of an effective, 
magnetic field-dependent number of channels. Predic- 
tions based on this approach are in good agreement with 
experiment. For multichannel samples with Af » 1 wc 
use the diagrammatic technique described in Refs. 38 and 
39. 

However, when interactions are present, this treatment 
should be modified. The approach which assumes that in 



■5 



a pointlike scatterer the interactions appear in the form of 
a self-consistent potential was introduced by Biittiker and 
co-authors on the basis of gauge-invariance and charge 
conservation. This (Hartree) approach neglects contribu- 
tions leading to Coulomb blockade (Fock terms), but is 
a good approximation for open systems. If the screening 
in the dot inside the medium with dielectric constant e is 
strong, r s = (fc^s) -1 = e 2 / (ehv-p) < 1, an RPA treat- 
ment of Coulomb interactions is sufficient. For large dots, 
L 3> a Bl the details of screening potential on the scale 
<~ a B arc not important and we can assign an electric 
potential U (r, t) defined by excess electrons at r, t at any 
point r of the sample. If additionally the number of bal- 
listic channels N is much smaller than the dimensionless 
conductance of a closed sample, gdot = i?Th/A » TV, the 
potential drops over the contacts and therefore in the 
interior of the dot it can be taken uniform ("zero-mode 
approximation"). 37 This potential shifts the bottom of 
the energy band in the dot and thus modifies the S- 
matrix. As a consequence, electrons with kinetic energy 
E have an electro-chemical potential E a — E—eV a in the 
contact a and E = E — elf in the dot. (We point out that 
we neglect the quantum pumping in the dot and conse- 
quently the 5-matrix depends only on one energy.) Re- 
cently, Brouwer, Lamacraft, and Flcnsbcrg demonstrated 
that this self-consistent approach gives the leading or- 
der in an expansion in the inverse number of channels 
1/N <C l. 41 Therefore, our analytical results present the 
leading order effect, valid for 1/N -C 1. 

In the self-consistent approach the influx of charge 
changes the internal electrical potential of the dot U (t) , 
which in turn affects the currents incoming through each 
conducting lead and/or redistributes charges among the 
nearby conductors (gates). Such capacitive coupling can 
often be estimated simply from the geometrical configu- 
ration. For example, the capacitance of a dot covered by 
a top gate at short distance d <C L is C ~ el? jd and a 
single quantum dot has C ~ sL. The ratio of charging 
energy E c ~ e 2 /C to mean level spacing A characterizes 
the interaction strength. It is proportional to the ratio 
of the smallest geometrical scale to the effective Bohr's 
radius, E c /A ~ min {d, L}/as- We refer to interactions 
as strong if E c » A and weak if E c <C A. 



IV. IMPORTANCE OF BIAS MODE 

We suppose for simplicity that at equilibrium the volt- 
ages V\ = Vi = Vo are set. In the following we consider 
the situation when the (single) gate voltage Vo is held 
fixed at its equilibrium value. Experiments can be per- 
formed in different bias modes, usually cither (i) with 
fixed drain voltage V2 or (ii) at fixed V\ + V2 (the vari- 
ations of the voltages at the contacts are equal in mag- 
nitude but opposite in sign). These different modes cor- 
respond to straight lines in the {Vi, V2} plane shown in 
Fig. 3. 

Let us consider the nonlinear current as a function 




FIG. 3: Depending on the bias mode, the experiment probes 
different transport properties. Plots present (left) linear and 
(right) nonlinear components of the current as functions of 
x = Vi — Vo and y = V2 — V0, the dashed curves correspond to 
equal currents. Thin line shows fixed V2 and I(V) is a function 
of source voltage V = Vi. Thick line corresponds to fixed 
Vi + V 2 , such that I(V) depends only on V = (Vi - V 2 )/a/2. 
Full and empty dots on the right figure correspond to the 
forward or reverse configurations shown in Fig. 1. 

I(x, y), where x = V\ — Vo and y = V2 — Vo are deviations 
of contacts voltages from equilibrium. For generality we 
consider below a situation when the linear combination 
— a;sin(f7 — 7r/4) +ycos(?7 — 7r/4) = is held fixed and the 
only variable is 

V = xcos(t7 — 7r/4) + y sm(r] — 7r/4). (3) 

This corresponds to a rotation of the original x, y axes 
such that the new coordinate axis V makes an angle 77 
with the y = —x line, as illustrated in Fig. 3. The 
value of 77 fully characterizes the bias mode. Now the 
two modes introduced above are simply (i) r] = 7r/4 which 
implies V = x; and (ii) ry = 0, which implies V = (x — 
y)/V% and corresponds to an asymmetric variation of the 
voltages. 

The linear current depends only on x — y (dashed lines 
on the left panel in Fig. 3 correspond to the lines of 
equal currents) and in any bias mode the measured linear 
current /n n is the same for a given x — y. If we consider 
the nonlinear current I as a function of x, y, it is by 
construction a bilinear function of x, y. As in the linear 
case the current must vanish if the voltages are the same 
and thus I = for x — y = 0. Therefore, the bilinear 
function must be of the form 

I = I [(x + y) cos <p+(x-y) sin 4>] (x - y) (4) 

with unknown (generally fluctuating) parameters Io and 
<fi G (— 7r/2, 7t/2]. It is important that the qualitative be- 
havior of I(x, y) depends on the interaction strength: one 
could expect that transport depends not only on voltages 
in the leads, but also on the internal noncquilibrium po- 
tential U of the sample. This potential can be found if 
potentials in all reservoirs and the nearby gate are known. 

In the limit of weak interactions the equilibrium point 
Vq is important, and if we reverse the bias voltage, 
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(V, 0) — > (0, V) the current is fully reversed, that is 
d 2 x I — —dyyl. For the current defined in Eq. (4) it is 
possible only when I oc (x — y) (x + y) => <j> = 0. Another 
way to see this is to use the usual expression for the total 
current in terms of scattering matrices. In this formula 
the current depends on the difference between Fermi dis- 
tributions in the leads oc f(e — ex) — f(e — ey), and its 
expansion up to the second order yields /"(e) (a; 2 — y 2 ). 
The lines of equal current are curved and directions 
1] = 0,±7r/2 correspond to zero current directions. Thus 
the dependence of current on the angle r] is strong. In 
addition this approach predicts that the current through 
a two-terminal sample is symmetric with respect to the 
magnetic flux inversion. 

In contrast, for strong interactions, the value of Vb is 
irrelevant and the nonequilibrium electrical potential U 
is independent of Vb . In this case current depends only on 
the voltage difference x — y and thus I cx (x — y) 2 => \<f>\ = 
it/2. The equal-current lines are straight and the picture 
is similar to the left plot in Fig. 3 for linear transport. 
Therefore we do not expect any nontrivial dependence of 
the nonlinear current on the choice of the bias mode. 

It is noteworthy that qualitative considerations can 
predict neither the sign, nor the magnitude of I®. The 
only general conclusion which we can make for a weakly 
or strongly interacting dot is I(x,y) cx x 2 — y 2 and 
I(x,y) cx (x — y) 2 , respectively. Experiments extract 
derivatives of I with respect to the applied voltages. Im- 
portantly, this derivative depends on the chosen direction 
77. The nonlinear current measured in this bias mode is 

I(rf) = IqV 2 cos?7sin(</> + rj). (5) 

The current is zero when 77 = —<f) and r\ = ±ir/2, and 
the bisectrix of the angle between the two zero-current 
directions at 77 = — </>/2+(7r/4)sgn maximizes d 2 I/dV 2 . 

Sometimes experiments extract information on non- 
linearity from measurements in different connections 
schematically shown in the central and right panels in 
Fig. 2: "Forward" connection corresponds to x = 
±V,y = 0, while "reverse" connection for the same 
voltage configuration corresponds to x = 0, y = ±V. 
The gate voltages V g are kept fixed. In Fig. 3 these 
forward and reverse points are indicated by black and 
white dots, respectively. To find the nonlinear conduc- 
tance Marlow et al. 8 and Lofgrcn et al. 19 determine the 
difference of conductance at these measurement points. 
Lofgrcn et al. 19,42 use the term "rigidity" for samples for 
which G f (V) = G r (-V) in the points /+ = (V, 0) and 
r~ = (0, —V). 42 Equation (4) gives the nonlinear contri- 
bution Q to the full conductance G/ jr (±V): 

Q cx I [(x - y) sin + (x + y) cos 0] . (6) 

Thus for a sample which is called rigid this implies 
Io cos 0^0. Since Iq = would mean that there is 
no second-order response, we must have cos <f> — * which 
is the case for samples with strong interaction. In other 
words, "rigidity" in samples which exhibit 0(V 2 ) current 
is equivalent to strong Coulomb interactions. 



On the other hand, comparison of data at another 
pair of points /+ = (V, 0) and r+ = (0, V) gives 
G f (V) — G r (V) cx Iq sin <j> and provides additional infor- 
mation about the two fluctuating quantities Iq, 4>- Refer- 
ence 19 expects that a Left-Right (LR)-symmetric system 
has Gf(V) = G r (V). Therefore rigid and LR-symmetric 
sample should necessarily have Jo — ► and thus could 
not exhibit a second-order current 0(V 2 ). This point is 
discussed more quantitatively in Sec. VA. 

It is important to note that to find the linear DC cur- 
rent one needs to know only x — y = Vi — V2 , while for 
the nonlinear current in general one needs two variables 
x = Vi — Vq, y = V2 — Vo or any independent pair of 
their linear combinations. The projection of the vector 
(Vi, V2, Vb) on the V\ + V2 + Vb =const plane uniquely 
defines the nonlinear current. This projection can be 
parametrized by the pair of Cartesian (x, y) or axial co- 
ordinates (V,r)). However, if in the experiment the volt- 
ages Vi j2 were fixed, this would not be enough to define 
(x, y) uniquely. In this case Ref. 19 points to the impor- 
tance of the reference point V . Indeed, one could arrive 
at the point with a given (Vi, V 2 ) from any equilibrium 
point and the measured current would depend on Vq. We 
prefer to characterize the measurement by the pair (V, r/) 
instead of three variables (Vi, V2, Vb) because of the sim- 
plicity of the final results. The weaker the interaction (or 
the stronger the capacitive coupling of the sample to the 
nearby gate) the more important the role of 77 chosen in 
experiment. 

We illustrate this important conclusion by quantitative 
results for nonlinear conductance Q cx d 2 I/dV 2 in the 
following sections. We point out that conductance with 
respect to the voltage difference V = V\ — V2 is often used, 
even when a linear combination V is actually varied in ex- 
periment. Voltages V and V are related, V = V/V2 COS J], 
and one can straightforwardly find d 2 I/dV 2 . 

V. GENERATION OF DC IN QUANTUM DOTS 

Now we quantify the qualitative arguments of Sec. IV 
and consider the more general situation of a DC current 
generated by an AC bias. If at first we neglect Coulomb 
interactions, the nonlinear DC current in response to 
the Fourier components Vp.^ — Vpe 1 ^ 13 of the AC voltages 
applied at the contacts (i = 1,...,M, can be expressed 
with the help of the DC-conductance matrix g a /3{e) of 
the dot at the energy e 28 

_ v s e 3 f fjs + M + f(e - M - 2/(g) 
a ~ h J {fuo) 2 

M 

x ^2g a p{e)\V fiiU \ 2 , (7) 
/3=i 

g a p(e) = ti[l a 5 a p-Si(e)TL a S(e)TL ]. (8) 

If we now include interactions using a self-consistent po- 
tential this formula is modified: 28 in Eq. (7) the 
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Fourier components of the voltages at all contacts are 
shifted down by the Fourier component of the internal 
potential — XJ W 



7 ^01 



G i g 7 (w) - ia;C s ' 



(9) 



v e 2 r 

G/3 7 (uj) = -j- J detv [l /3 l 7 - l 7 5 t (e)1^5(e + /^)] 



/(g)-/(e + M i 



(10) 



In Eq. (9) the index 7 runs not only over real leads 
1,...,M, but also over all gates gi. However, when 
7 G {ff?} the AC conductance G/3 7 (w) is absent and only 
capacitive coupling iwC 7 remains in the numerator. We 
point out that the matrix G{u>) of dynamical AC con- 
ductance at frequency ui given in Eq. (10) should not 
be confused with the degenerate matrix g(e) of energy- 
dependent DC conductances of electrons with kinetic en- 
ergy e given in Eq. (8). 

The results of Ref. 28 can be expressed in terms of the 
DC conductances g a p and frequency-dependent charac- 
teristic potentials u 7 , 

^e 3 f /(e + M + /(£-M-2/(g) 

' /< J ' ~ (M 2 

x ^2g af3 (e)Re u 7 |Va, w - V JtbJ \ 2 . (11) 

0i 

Here Re u 7 stands for the real part of u~ n which is in 
general a complex quantity. In contrast to Eq. (7), Eq. 
(11) is expressed via differences of voltages applied to 
all present conductors. Therefore, the current is gauge- 
invariant. The charge conservation, ^ Q I a = 0, is obvi- 
ous from Eq. (8). 

From this point on we consider Eq. (11), a specific 
expression of Eq. (1), in detail for several regimes. In 
Sec. VA we discuss the nonlinear current due to DC 
applied voltages (previously considered in Ref. 43) and 
the importance of different bias modes in experiments in 
two-terminal quantum dots. In Sec. V B we consider the 
frequency dependence of Q s {uj) and G a {w). 



A. Nonlinearity in quantum dots 

In the static limit 43 fojj/T — > the integrand in the first 
line of Eq. (11) simplifies to f"(e) and for Hlu/NA -> 
the derivatives u 7 are real and expressed via subtraces of 
the Hermitian Wigner-Smith matrix d E S / (2Tri) 44 ' 45 



I n = 



-v s e" 



£ [ f'{e)deg' ap {e)u 1 {V - V^f, (12) 

= C 7 Me 2 - J def (g)tr l^d e S/(2m) 
Ul Cs/^e 2 - J def'(s)tvSW s S/(2ni) ' 1 ' 



For a two-terminal sample the nonlinear current through 
the first lead is 



h = 



f'(s)g' n (e)de 



iV 2 -V gl f] + (u 2 -u 1 ){V 1 -V 2 f 



(14) 



The characteristic potentials in the last term of Eq. (14) 
are sensitive to the asymmetry of the contacts. Indeed, 
in a strongly interacting dot u g i = and u 2 — u\ w (N 2 — 
Ni)/N. The current magnitude grows with asymmetry 
due to the last term in Eq. (14). On the other hand, the 
sign of I\ is random because of quantum fluctuations of 
g'n around zero. 46 As a consequence, if in an experiment 
the Fermi level is shifted by S/ip ~ NA/2ir (or the shape 
of the dot is changed) the sign of nonlinearity can be 
inverted. 

Different modes of bias having been discussed in Sec. 
IV, we concentrate here on the (anti)symmetric conduc- 
tances through the quantum dot at fixed gate voltages. 
When the reservoir voltages are varied in the 77 direction, 
the nonlinear current is given by the expression 



I = 



-2v,e 3 



//'(eKiOOcfe [(!-«! - w 2 )sin?7 



+ {u 2 — ui) cos 77] cos r]V 2 



(15) 



and one can define exactly the unknown parameters Iq, 4> 
which we introduced in the qualitative argument leading 
to Eq. (5). Depending on 77 one measures different linear 
combinations of conductances. If we consider conduc- 
tances d 2 I/2dV 2 in units of v s e z jh, Eqs. (2) and (15) 
yield 

Q 2tt cos 2 77 J dedef{e)f{e)xi (e)X2,a(s) (£) 

ya - s ~ A 2 [C s /(eV s ) - J def (e)trSW £ S/(2m)} 1 ' 

expressed in terms of fluctuating functions \ an d a trace- 
less matrix A = (N 2 /N)\i - (N 1 /N)l 2 : 



Xi (e) = (A/2^)a £ trA5 t A5, 
X2,a(e) = (*A/27r)trA[5t,a e 5], 



(17) 
(18) 



X2, s (e) = A ( - + 



N 



2m 



(19) 



Standard calculations using the Wigner-Smith and/or S- 
matrix averaging 38 ' 39 ' 47 yield (Q a ) = (Q s ) = 0. This re- 
sult signifies that the nonlinear current through a quan- 
tum dot is indeed a quantum effect. As a consequence 
the size of the measured nonlinearity must be evaluated 
from correlations of G a ,G s - 

The functions Xi( £ >$) an( i X2,a/s( £r '> arc uncorre- 
cted, and their autocorrelations 16 readily allow one to 
find statistical properties of G a ,s- Our results can be ex- 
pressed in terms of diffuson T> or coopcron C in a time 
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representation, exp(— t/tjj) and exp(— t/tq)- Both can 
be introduced using the 5-matrix correlators 48 (correla- 
tions of retarded and advanced Green functions lead to 
the same expression up to a normalization constant 37 ). 
We have 



ier /h 



and actually would be best for an accurate measurement 
of linear transport. Indeed, this regime is not affected 
by the fluctuations of capacitive coupling uo of the dot 
with the nearby gate and thus minimizes fluctuations of 
G s around 0. 

Fluctuations of Ga,G s are given by different expres- 
sions, see the first line of Eq. (23), where the first term 



(Sy(T, $')) - {e-^SikSjt+e-^SuSjk) is duc to (X 2 ,a) ° r ^J- Importantly, (x 2 , s ) contains 



(20) 



1~C,T> = 



N c 

N c ,vA' \Nt> 



4$2 



L 2 A' 



We also introduce the electrochemical capacitance C M 49 
which relates the non-quantized mesoscopically averaged 
excess charge (Q) in the dot in response to small shift 
of the voltages 5V at all gates. In addition the charge 
relaxation time trc of the dot is conveniently introduced 
by this electrochemical capacitance and the total contact 
resistance, 



C„. 



(SQ) 
5V 



l + C s A/(^e 2 ) : 



trc = 



v s Xe 2 ' 



(22) 



The denominator of Eq. (16) is a self-averaging quantity, 
((...) 2 ) = ((...)) 2 = A 2 (C S /C M ) 2 . Using the diffusons 
and cooperons defined in Eq. (20) we find the following 
correlations of Go, and Gs- 



(&(*)&(*')> 
x hcos 2 ,-^ 



N*N* 



AT 



dr. 



X = 



\2h 2 ) J sinh 2 7rTr/ft' 
N 2 /C tan77 + C 2 -Ci 



2N X N 2 



,e 2 /A 



+ 



N 2 - Ni 
N 



(23) 
(24) 
(25) 



There are two very different contributions to Eq. (23), 
Tc,t> due to quantum interference and X defined by the 
classical response of the internal potential to external 
voltage. The terms denoted by Tc,t> are sensitive to 
temperature, magnetic field, and decoherence. Asymp- 
totical values of T in the low temperature, T <C H/t\, or 
high temperature limits, T » fo/r\, are T\ — ► 1/N 2 — 
(r x A/h) 2 and T x -» A/(12TA A ) = T X A 2 /(12hT), re- 
spectively. 

The term denoted by X and given by Eq. (25) contains 
only quantities specifying the geometry of the sample and 
gates and the bias mode. In a real experiment the cou- 
pling due to capacitances d j2 is usually stronger then 
that of the external gates, Ci l2 3> Co- Symmctrization of 
the circuit C\ = Ci can diminish the value of X. If in ad- 
dition Ni — N 2 and 77 = (used in the experiments 9 ' 10 ) 
we have X — > 0. Thus such a symmetric setup and bias 
mode minimize the fluctuations of the nonlinear current 



both quantum T\ < 1/A^ and classical X contributions. 
If the classical term dominates, X ^> 1/N 2 , the current 
is mostly symmetric, (? 2 » Q\. This could be expected 
either for a weakly interacting dot or a very asymmetric 
setup, N\ 7^ 7V2 . 16 However, if the classical term is re- 
duced due to, e.g., the bias mode, the fluctuations of Ga 
and Gs become comparable. This experimentally impor- 
tant conclusion remains valid for any interaction strength. 
(Particularly, it leads to a very wide distribution of the 
Aharonov-Bohm phase considered in Sec. VI.) 

Experiments of Zumbiihl et al. 10 and Leturcq et al. 9 
are performed in this regime when rj = and X — > 0. 
Data in Ref. 10 demonstrate that the part of the total 
current symmetrized with respect to magnetic field is by 
far dominated by linear conductance. From Eq. (23) 
we expect mesoscopic fluctuations in linear conductance 
to be ~ iV 2 times larger then those of (? S A. Thus only 
when the number of channels is decreased will the non- 
linear Gs become noticeable. A clear observation of Gs 
without linear transport contribution was performed in 
a DC Aharonov-Bohm experiment by Angers et al. 11 in 
the mode n = ±7r/4 (only one contact voltage was var- 
ied). This allowed to evaluate the interaction strength 
from the ratio of GsjGa- 

Experiments of Marlow et al. 8 and Lofgren et al. 19 
measure the full two-terminal conductance and extract 
nonlinear conductance properties related to various spa- 
tial symmetries of the dot. Although the current through 
a weakly interacting sample is field-symmetric, this is not 
true in general. Samples of Ref. 19 differ in "rigidity" and 
degree of symmetry. Rigid samples, u — » 0, with Left- 
Right (LR) and Up-Down (UD)-symmctry should have 
(u2 — «i) s = and (u 2 — u\) a = respectively, according 
to the expectations of Lofgren et al. 19 (indices s and a 
mean the symmetric and antisymmetric part in magnetic 
field). 

Duc to quantum fluctuations, in experiment none of 
these symmetry-relations can be exactly fulfilled, see 
Eq. (23). According to Eq. (15), the difference in 
the full conductances g — {h/v s e 2 )I /V measured be- 
tween different points probes different characteristic po- 
tentials. Reference 19 defines three differences 5i,ii,iii for 
three pairs of points in the forward and reverse con- 
nection discussed after Eq. (5). Using Eq. (23) we 
find (i) gi = gf(V,B) - g r (-V,B) cx u , (ii) g i{ = 
9f( v , B ) ~ 9f(V,-B) cx (u 2 - ui) , and (iii) g m = 
gf(V, B) — gf(—V, —B) oc (u + u 2 —Ui) s . The ensemble 
average of these differences vanishes and their fluctua- 
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tions for Ci i2 — 0, N\ = N 2 are given by 









H 


it J 






neV 
2A~~C 



where X = 2(C/C* M - l) 2 is found from Eq. (25) at 
i] = ±7r/4. In weakly interacting dots C^/C — > and 
only magnetic-field symmetric signals (jj and t/iii survive. 
In strongly interacting ("rigid") dots C^/C — > 1 and <7;i 
becomes similar to gm. We point out that even if the 
rigid samples are made symmetric with respect to Left- 
Right inversion, the quantum fluctuations of the sample 
properties are unavoidable and gf { ^ at $ ^ 0. For 
high magnetic fields and arbitrary interactions Tc —> 
and experiment should observe gf + gf = gf Vl . Clearly 
fluctuations exist also for large magnetic fields beyond 
the range of applicabilty of RMT. Experimental data (see 
inset of Fig. 6 in Ref. 19) show that gf + gf ~ gf-. It 
is hard to make a quantitative comparison with Rcfs. 19 
and 8, since the quantum fluctuations in the nonlinear 
conductance exist possibly on the background of classical 
effects due to macroscopic symmetries. We expect that 
quantum effects become more pronounced as contacts are 
narrowed. 

To conclude this subsection we briefly discuss here the 
case of a macroscopically asymmetric setup. If the ex- 
periment were aimed to measure large Q a compared to 
Q Sl one would try to minimize Q s by adjusting the setup. 
Such a procedure minimizes the value of X in Eq. (25). 
For Ci ; 2 = 0,77 — 7r /4 the role of asymmetric contacts 
N\ 7^ N 2 was discussed in Ref. 16. Analogously, one 
could consider a more general case of C\_ 2 and an ar- 
bitrary bias mode 77. This is especially important if the 
difference C\ 7^ Ci can not be neglected due to occasional 
loss of contact symmetry. 

The results of an experiment could also be affected 
by the presence of classical resistance loads 7*1,2 be- 
tween macroscopic reservoirs and the dot (shown in 
Fig. 2). Swapping of such resistances in the ex- 
periment, when connection is switched between "for- 
ward" and "reverse" 19 affects the voltage division be- 
tween loads. If we assume the capacitive connection of 
the dot and reservoirs is still the same, the modification 
of the expression for u 7 in Eq. (9) is straightforward, 
Y,p G Pi -> E/3 G Hil ( 1 + r i E/3 G Pi)- Naturally, at large 
r 7 , (2e 2 A 7 //i)r 7 3> 1, the main drop of the voltage oc- 
curs over the resistor r 7 and not over the QPCs. As a 
consequence, if n,2 7^ 0, values of wi j2 can become un- 
equal due to n 7^ r 2 and this leads to the classical circuit 
asymmetry which we do not consider here. 



B. Rectification in quantum dots 

Here we consider the DC generated by a quantum dot 
subject to an AC bias at the frequency to. In experi- 
ment at high bias frequency lut^ > 1 current is usually 



measured at zero frequency. In contrast, at small bias 
frequency lut^ <C 1 higher harmonics (for instance the 
second harmonic 2lu) can be measured. However, up to 
corrections small due to lut^ <C 1, the second harmonic 
is just equal to the rectified current, / 2w ~ h- There- 
fore, to leading order, our results for the rectified current 
describe both experiments. 

Generally, there are several important time-scales char- 
acteristic for time-dependent problems in chaotic quan- 
tum dots. To see how they appear let us first consider 
frequency-dependent linear transport of noninteracting 
electrons. Its statistics usually depend only on the flux- 
dependent time scales tq,v, see Eq. (21). If we consider 
an analog of UCF (G 2 ($)) for the frequency-dependent 
conductances introduced in Eq. (10), we find 



v s e 2 Nj_N 2 
h N 

r -l)(l- 



E 

\=c,v 



( AT\ 2 f T X dr e- T l T *(e lL 
\2^J J lulu' [1 - i(uj + uj')t x / 2} sinh 2 ttTt/H 



• (26) 



The presence of iu)T\ in the diffuson and cooperon contri- 
butions in the second line of Eq. (26) is due to the energy 
dependence of the scattering matrix S(e), which usually 
brings up imaginary corrections to the matrix-element 
correlators. 

In a DC-problem lu — > it is usually useful to intro- 
duce a dimensionless number of channels Nc,t> modified 
by the magnetic field, see Eq. (21). In this limit at 
T — > the integration in Eq. (26) becomes straightfor- 
ward and summation is then performed over A^ 2 . For 
equal magnetic fields, $ = we have N-p — N, but 
Nc is strongly modified by large fields, Nc — > 00, which 
suppresses the weak localization correction and dimin- 
ishes UCF. However, for an AC problem (especially for 
LUTd > 1) it is more convenient to express results in terms 
of dimensionless quantities lutc,lutx>. For example, from 
Eq. (26) the statistics of conductance can be easily eval- 
uated: (|G(w,<I>)| 2 )/(G(0,$) 2 ) ~ l/(wT d ) and the real 
and imaginary parts of conductance are similar and un- 
correlated at high frequency lut^ 3> 1. 

Inclusion of interactions introduces an (additional) de- 
pendence on trCj the charge-relaxation time defined in 
Eq. (22). To leading order in 1/N < 1 the effect of 
interactions is often to substitute Td — > trc in the non- 
interacting results, e.g., for the linear conductance 49 or 
shot noise. 45 ' 50 Interestingly, the subleading corrections 
depend on both trc and Td, e.g., in the weak localization 
correction in the absence of magnetic field. 49 When the 
magnetic field is increased to values which finally break 
time-reversal symmetry, the appearance of different time 
scales rc,x> is expected, see e.g., Eq. (26). Therefore at 
intermediate magnetic fields, when td ^ t~c and the in- 
teractions taken into account, trc 7^ T d, the solution of 
an AC problem is expected to show a complicated depen- 
dence on all these time scales. 

Indeed, if we consider the rectified current such an 
interplay between tc,t> and trc does appear. We find 



10 



(Ga) = (Gs) =0 and present below results for correla- 
tions of G a and Gs'- 



(&($)&(*')> 



X^G,C(U)+^(U)] 



Fu,t>{w) - Fu,c{u) 
47rcos 2 ry iVfiVf 



N 6 



.(27) 



Here the functions Fu{uj), Tg{<~°) are finite-frequency 
generalizations of Eq. (24) 



(try j 



T\dr 



e -T/T X gj n 2 UT j 2 1 



\h 2 u) J ' A "" 2sinh 2 7rTr/ft 1 + 
1 + !wr RC e iwr 



x 1 + Re 



AT 



ZWTrc 1 - ibJT\ / 
2dr e- T / T >- sin 4 wr/2 
r A sinh 2 7rTT/n 



(28) 



(29) 



The subscripts U(G) of .F^g) (w) illustrate the origin of 
these functions: they result from averaging of different 
scattering properties over the energy band defined by 
max {hu,T,h/Tc(v)}- The function J 7 u(^>) is a charac- 
teristic of the internal potential U u , see Eq. (9). The 
function ^Fq(u>) results from the energy averaging of the 
DC conductance g(e). Such averaging appears because 
both G(ui) defined in Eq. (10) and g(e) in Eq. (11) are 
coupled to the Fermi distribution. 
The function X(u) is 



N 2 



2N 1 N 2 
N 2 



(C tan 77 + C 2 -C 1 )(l + w 2 r d T RC ) 



(l+c. 2 r 2 c )^e 2 /A 



+ 



7V(1+u, 2 t 2 c ) 



(30) 



and in the static limit u — ► it is given by Eq. (25). We 
point out that when the interactions are negligible, E c ~ 
e 2 /C <C A, the role of the bias mode is significant. A 
quantum dot with (fully broken) time-reversal symmetry 
can be labeled by Dyson symmetry parameter ((3 = 2) 
(3=1. When the setup is ideal, C\, 2 — 0, and r\ ^ 0, the 
fluctuations of Q a , Gs at large frequencies uot^ >• 1 are 



SGs 

SGa 



(Gt) 



2\l/2 



N1N2 

N 2 



2 7T \ 1/2 2sin277 

(3 LUTd J Hll> 



iViiV 2 \ 3/2 v s e 2 cos 2 7? 

N 2 ) 2C H 2 LU^T d 



(3 = 2. 



(31) 
(32) 



In chaotic quantum dots the role of the Thouless energy 
-Erh of the open systems is often played by the escape rate 
h/ra- If we take this into account, our result (31) qual- 
itatively agrees with that of Falko and Khmelnitskii 20 
obtained for open diffusive metallic junctions. However, 
when 77 — > 0, the fluctuations of G s are much smaller and 
for I TV"! — 7V 2 1 <C they become comparable with those 
of the antisymmetric conductance (32). 



TABLE I: Asymptotes of Fu,\(u), X(w) at T -» 





Adiabatic 


Intermediate 


High 


Function 


co « r- 1 


r- 1 « w « r R( ! 


WTRC > 1 


Fu(w) x JV| 


1 


7r/(4wr A ) 


tt/^tat^c) 




1 


(tt/Wa) 3 


X{u) 


2 tan 2 r?(r R c/r d ) 2 


2 tan 2 T) 



However, very often experiments are performed in sam- 
ples, where the interaction is not weak. Since Ga and Gs 
are comparable for strong Coulomb interactions in the 
static limit ui — > 0, 16 we concentrate here on this exper- 
imentally relevant regime of A/E c <~ TRc/Va <C 1 and 
take an ideal symmetric setup, N\ — N 2 and Ci,2 = 0. 
Then we have 



X(to) « 2 tan^ 77 — 

W c +w 2 t rc 



(33) 



Below we consider in detail the case r? 7^ and how this 
bias mode affects the behavior of £ 2 (w). Several fre- 
quency regimes can be separated: adiabatic lot\ <C 1, 
intermediate, where 1/t\ <C lj <C 1/trc, and high fre- 
quencies wtrc ^ 1- The asymptotes of the functions 
defined in Eqs. (28), (29), and (33) in these regimes are 
presented in Tabic I for reference. 

For adiabatic frequencies lot\ <C 1 the integrands in 
Eqs. (28) and (29) do not oscillate on the short time 
scale t\. At such small frequencies Tu(w) = Tg{oj) are 
equal to J 7 of Eq. (24) and J(w) cx (r RC /Td) 2 <C 1 can be 
neglected. This is essentially the zero frequency regime 
considered before for nonlinear DC transport. 

As the frequency grows, an intermediate regime is 
reached when max {T,h/r\} <C fou <C /j/trc and 
!Fu(io), Tg{u) start to differ. The scattering properties 
at large energy difference ftw 3> h/r\ arc uncorrelated 
and the response of the dot is randomized. Therefore 
both the conductance averaged over a large energy win- 
dow hu> and the response of the internal potential U u to 
the AC voltage at lot^ >• 1 are strongly suppressed, see 
Table I. As a result, if X(lj) is still negligible, both Ga 
and Gs decrease with growing frequency as 1/w 4 . 

One could expect that interactions qualitatively 
change the behavior of G a ,G s when the frequencies be- 
come comparable to 1/trc ~ NE c /h, the scale defined 
by the interaction strength. At such frequencies the re- 
sponse of a dot to the potentials at the contacts is not 
resistive as occurs at low frequencies, but mostly ca- 
pacitive. If the frequency is high, wtrc > 1, we have 
Re Mi,2 — > and the function Tjj in Eq. (28) is sup- 
pressed ~ 1/(wtrc) 2 - As a consequence, Q\ is suppressed 
stronger then 1 /lj 4 and goes as 1/ui 6 at cjtrc > 1 . How- 
ever, a more important signature of this capacitive cou- 
pling is the growth of X(lj) in Eq. (33), which affects 
Gl 

To see the role of this growth we consider now suffi- 
ciently large fields <& = <!>' when only the diffuson contri- 
bution survives. The growth of X(u>) in Eq. (33) reflects 
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TABLE II: Estimates for (gl) and (5?), {z = u)t a ) 




Function 


Adiabatic 
2 < 1 


Intermediate 

1 < z < r d /r R c 


High 

z > r d /r R c 




1 


2- 4 


^2- 6 


T TtC 


1 


(1 + ^)2,-1 


4-z-* 



we predict a monotonic decrease of (G1{oj)), this growth 
could be the result of quantum pumping or an interfer- 
ence of the pumping and rectification (both effects were 
neglected here). 



FIG. 4: Zero-temperature large-field fluctuations of Q a (w) 
(dashed) and Q s (u) (solid curve) in units of (tv/4AN 2 ) 2 for 
the bias mode r\ = 7r/4. Data are presented in the log-log scale 
at iVi.,2 = 5 and Tnc/Td = 0.05. The asymptotes Q\ oc lo~ 6 
and C/s oc oj~ 3 are different due to r? ^ 0, see Eqs. (34) and 
(35). 



enhanced sensitivity of the internal potential to the 
gate voltage, X(u) oc (tanr/Re u ) 2 . At high frequencies 
the impedance of the capacitor C becomes negligible and 
therefore the internal potential follows the gate voltage 
and not the reservoir voltages, uq — ► l,iti,2 — > 0. En- 
hanced from its small static value trc /ra to 1 at large fre- 
quencies, such coupling affects the fluctuations of Q s (oj) 
if 7] ^ 0. The situation is somewhat similar to the weak 
interaction limit, when the coupling with nearby gates 
was strong, u — > l,«i,2 -C 1, and lead to Q s >• Q a . 

The fluctuations of Q a (uj),Q s (uj) for ut^ 3> 1 can be 
evaluated: 



(M 4 (1+U 2 T, 2 



(34) 



RCV 



r 2/ \ r 2, \ , (r RC tan?7[l +u; 2 T d T R c]) 2 . . 
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Fluctuations of <? 2 (w) and <? 2 (w) demonstrate qualita- 
tively different behavior, which we illustrate in Fig. 4. 
Indeed, at sufficiently high frequencies, the dependence of 
X(ui) on uj makes the last term in Eq. (35) dominant. At 
wt rc 3> 1 the asymptotes of G 2 a oc 1/ui 6 and Gl oc 1/oj 3 
become different due to the presence of the second term 
in Eq. (35). These results show that for nonadiabatic fre- 
quencies of the external bias the DC current strongly de- 
pends on the bias mode r\. We predict that the magnetic 
field asymmetry of the rectified current, noticeable at 
small frequencies, might become suppressed for large fre- 
quencies, when the symmetrized component dominates 
due to the presence of capacitive coupling. For conve- 
nience, the low-temperature estimates for ((J 2 ) and (Gl) 
for r\ ^ 0, <f> > $ c are collected in Table II. 

It is noteworthy that a recent experiment in AB rings 12 
finds that Q(oj,$ = 0) grows with frequency until u> <~ 
2i?xh and then decreases <~ l/a> 3//2 ,aj — > oo. While 



VI. PHASE OF AHARONOV-BOHM 
OSCILLATIONS 

In this section we consider nonlinear transport through 
a chaotic Aharonov-Bohm (AB) ring. The nonlinear con- 
ductance Q exhibits periodic AB oscillations and non- 
periodic fluctuations, similarly to the linear conductance 
G. However, since Coulomb interactions produce asym- 
metry of G with respect to magnetic field inversion, the 
phase of these oscillations is not pinned to (mod)7r. As 
a quantum effect this AB phase is characterized by a 
mesoscopic distribution. The width of this distribution 
represents a typical fluctuation. We first discuss what 
kind of distribution could be expected in a chaotic AB 
ring and then calculate the fluctuation of the AB phase. 

Let us assume that Q as a function of magnetic flux 
$ can be expanded into the series of well-defined Fourier 
harmonics similarly to the linear conductance G: 
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• (36) 



The phase 5 of the main (first) harmonic <&o = hc/e is 
obtained from the ratio of the (anti) symmetrized con- 
ductances defined in Eq. (2) 



tan S 



fd$ CXp(27ri$/$ )^a($) 

J d$ cxp(27ra$/$ )£ s ($)' 



(37) 



We can not find the full mesoscopic distribution of the 
phase P{6). We can gain some insight in the behavior 
of this phase by investigating a similar quantity, namely, 
the asymmetry parameter A = Ga/Gs considered previ- 
ously for chaotic dots. 16 Based on Eq. (37) we argue that 
the statistical properties of arctan*4 and the AB phase 6 
should be similar. 

In quantum dots the parameter A is given by the ratio 
A = Ga/Gs = X2a/X2s, see Eqs. (16), (18), and (19). 
The functions X2a,2s at T ^= are convolved separately 
with /'(e), and at T — (which we consider below) they 
are evaluated at the Fermi energy. The properties of 
Xia.2s and the dependence of \2s on the bias mode were 
described after Eq. (25). The function \2s can have 
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FIG. 5: Mesoscopic distribution P{(j>) of (j> — arctan Q a /Qs- 
(Main plot) If the contacts are asymmetric (bold curve, N = 
16, Nl = 4) the distribution is narrow, while for symmetric 
contacts (dashed, N = 16, Nl = 8) it is almost uniform. As 
shown in the inset, for symmetric contacts at large N the 
distribution becomes uniform, compare bold curve for N = 2 
and dashed for N = 24. 

a nonzero (classical) average (X2s) ~ X 1 ' 2 defined by 
the interaction strength, geometry of the setup, and the 
bias mode r/. Since (x2,a) = and the fluctuations of 
X2a,2s are small as 1/N 2 , the mesoscopic distribution of 
arctan A is narrow and concentrated close to 0. How- 
ever, (x2s) = is possible if X — > 0, e.g., for symmetric 
contacts and the bias mode r) = 0. In this case, the 
distribution of arctan A becomes wide regardless of the 
interaction strength. 

The role of the classical contribution on the shape of 
P(arctan*4) is demonstrated in the main plot in Fig. 5 
for 7] = 0, where the distributions for asymmetric, N L = 
4, N = 16, and symmetric contacts, N L — 8, N = 16, 
are presented. While the distribution is almost uniform, 
when the classical contribution X is absent, it is highly 
peaked near zero when X dominates. If X is absent, the 
correlations between Q a and Q s are significant at small N. 
This leads to a nonuniform distribution of P(arctan„4), 
which is peaked at and 7r/2 when N = 2, see the in- 
set in Fig. 5. When A" grows, the correlations between 
Q a and Q s vanish and therefore the distribution becomes 
uniform. Such a distribution could be easily obtained if 
we make the natural assumption that Q a , Q s are indepen- 
dent and distributed by the Gaussian law with the same 
width. 

These numerics were performed for r\ — 0, when the 
mesoscopic distribution of A = G a /Gs becomes insensi- 
tive to the interaction strength. The role of interactions 
appears only if rj ^ 0, when the classical contribution X 
becomes dominant. Similarly, we expect that the distri- 
bution of the phase of AB oscillations is also strongly af- 
fected by the bias mode. If the bias mode is chosen such 
that the classical contribution X vanishes, the phase 5 
strongly fluctuates even for weak interactions. It would 
be very interesting to check this surprising conclusion ex- 




FIG. 6: Model of a chaotic Aharonov-Bohm ring with N = 
Ni + N2 channels. The model consists of a quantum dot with 
M channels combined with a ballistic arm with A3 = A4 = 
(M - N) /2 channels. 



perimentally. 

Let us now consider the fluctuations of the AB phase. 
Since the scattering theory turned out to be very use- 
ful for the discussion of the nonlinear/ rectified current 
through a chaotic quantum dot, we extend this theory to 
rings. We make two key assumptions (discussed in the 
Appendix in more detail) that the magnetic flux through 
the annulus of the ring is smaller then the flux quantum 
$0 and that the mean free path I, the radius R, the width 
of the ring W and the contacts W c satisfy the condition 
tt 2 IW > 2RW C . In this case the RMT can be applied 
to such chaotic rings as well. Unlike the experiments on 
large open rings with high aspect ratio R/W ^> l ; 2 ' 22 ~ 25 
the recent experiments 9,11,12 arc performed in rings of 
submicron size, which are effectively zero-dimensional. 
The treatment of such rings is similar to chaotic quan- 
tum dots, and the fluctuations of Q s , Q a can be expressed 
in terms of the diffuson T> and the cooperon C, see Eq. 
(23). The only problem is to find the expression for the 
effective number of channels as a function of magnetic 
field, similar to Eq. (21). 

The model we propose for a chaotic AB ring combines 
chaos and a ring geometry: a chaotic dot is attached to 
a long ballistic arm which serves to include an AB flux 
large compared to the fraction of the flux through the 
sample. This model is shown in Fig. 6, where the ring 
with A^ = A^i + A^ 2 ballistic channels in the contacts 1,2 
is modeled by a dot with M > N channels and a ballistic 
arm with A3 = N4 — (M — N)/2 channels in contacts 
3,4. The parameter p=l- N/M, the ratio of A3 + N 4 
to the total number of channels M, can vary between 
when the arm is much narrower then the contacts and 1 
in the opposite limit. The electronic phase is randomized 
in the quantum dot, but when electrons propagate along 
the arm their phase is determined by the geometry and 
applied magnetic field. This model is a reasonable ap- 
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proximation for the real experiment, it takes into account 
the long time spent by electron inside the ring and the 
randomness of its motion. The discussion of the model 
and the details of calculation of C,T> are presented in the 
Appendix. 

In experiment the Fourier transform is often taken over 
the total flux (or applied magnetic field) and the flux 
through the hole $/j cannot be separated from the flux 
through the dot Then the dependence of the diffuson 
and cooperon on magnetic field is non-periodic, which is 
indeed observed in the form of nonperiodic fluctuations 
in the (non-)linear conductance and phase slips of AB 
oscillations. A possible weakness of this model is in its 
spatial separation of chaotic scattering and the main part 
of magnetic field, but in the limit when the arm is much 
wider then the contacts 1 and 2 such a separation is not 
important and the averaged properties of AB oscillation 
phase become independent of the arm's width. 

If the flux $d through the dot is much smaller then 
the flux <&h through the hole, the nonperiodic fluctua- 
tions and the periodic AB oscillations are well-separated, 
which is usually the case in experiment. 9 ' 11 In view of this 
separation we can neglect the flux through the chaotic 
dot, &d <C to find the statistics of the AB phase. 
We assume that the averaging is taken over a magnetic 
field range containing many AB oscillations but still small 
compared to the characteristic field of the nonperiodic 
fluctuations. In such a simplified model of a chaotic AB 
ring Nc,t> are given by Eq. (A6) with <f> d = 0, $/j = <f> 
and the parameter p = (M — N)/M. The effective num- 
ber of channels is 



{£}-"('->-^'i- <»> 



Using this expression for the Nc,v we can evaluate the 
quantum fluctuations of linear conductance in AB rings. 
At low temperature a typical fluctuation of G at $ = is 
SG = ^{N^/N 2 )^ 2 /h) and the amplitude of AB 
oscillations is 8G\ B ~ SG 2 p(l - p) 1/2 /(l + p) 3/2 , which 
reaches maximum when the widths of the arm and the 
contacts are equal, p = 1/2. 

The first two moments of tan S can be found analyt- 
ically. It is zero on average, (tan 5) = 0, since the nu- 
merator and denominator in Eq. (37) are independent 
random quantities. Equation (38) show that T>,C are the 
same functions of <& ± <&', so that all necessary ingredi- 
ents can be expressed in terms of the functions Tu,v and 
T G .v of $ - $' defined in Eqs. (28) and (29). We omit 
now the index T> for brevity and denote the average over 
magnetic field by (...) = J*° (...)d<f>/$ , to find 



(tan 2 8) = T v (§)T G (*) cos $ + T G ■ T v ($) cos $ 



-Tu -:F G ($)cos$ / (Tu{§) + A(w))J r G ($)cos$ 
+ ^4-^g(*) cos $ + FuW) cos * . ( 39 ) 



where the function X(lu) is defined by the setup geometry 
in Eq. (33). Again, in the static limit ui — > we have 
Tu = T G =T and X defined by Eqs. (24) and (25). In 
this case Eq. (39) can be rewritten as 
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In the limits of high, T 3> NA/2n, and low temperature, 
T <C NA/2-k the asymptotical values of (tan 2 8} are 
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Very important is the case of symmetric contacts, N\ — 
N2, and antisymmetric bias mode, 77 = 0, which is used in 
Ref. 9. Then X vanishes and the average tan 2 5 becomes 
independent of interaction strength and as a function of T 
it is very weak. That is not the case if r\ ^ 0, for example, 
when only one of the voltages changes, r\ — zkir/4. 11 ' 12 
Then the statistics of the AB phase becomes temperature 
and interaction dependent due to the presence of X. 

The limit M ^> N corresponds to a uniformly chaotic 
ring, which we suppose to be closer to the experimental 
situation. Then the dependence on M drops out and the 
high/low temperature asymptotic read 



1 



(tan 2 S) 



X (3NT/A, T»NA/2,, 

\ N 2 /5, T -c NA/2n V ' 



This result clearly demonstrates that the phase of the 
oscillations is expected to deviate strongly from 0, es- 
pecially if the temperature is low and the number 
of channels in the contacts is diminished. The tem- 
perature is taken into account only in the form of 
temperature-averaging and the dephasing (previously 
considered for nonlinear transport of noninteracting elec- 
trons in Refs. 51-53) is not included. 

We expect our model for chaotic AB rings to work both 
for experiments at small frequencies 9 ' 11 and for large 
frequencies. 12 Similarly to quantum dots, the generaliza- 
tion on the finite-frequency case is obvious, if we use Eq. 
(33). Even in cases where RMT cannot be assured to be 
valid for open diffusive rings, the dependence of the AB 
phase on interaction strength, temperature, and number 
of external channels given by Eq. (41) should be correct 
qualitatively. 

The experiment of Leturcq et al. 9 is performed in a 
bias mode r/ = when X = 0. Then Eq. (41) gives 
(tan 2 S) = 1. The phase of the oscillations is evaluated 
from data according to Eq. (37) over a large range of 
fields. In experiment the AB phase is varied continu- 
ously as a function of the gate voltage at one of the arms 
of the ring. The data demonstrate that the phase 5 in- 
deed changes in a wide range and is usually far from 
0. This substantiates our conclusion that in the mode 
when the classical contribution is minimized, X — > 0, the 
mesoscopic distribution of 8 is very wide. 
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Experiment of Angers et al. 11 varies voltage in a dif- 
ferent way, T] = n/ 4, and therefore has 1^0. We would 
expect the phase 5(mod)7r to take values closer to and 
the antisymmetric component of the oscillations be rela- 
tively smaller even for large fields. Although phases close 
to are indeed observed, the field averaging is taken only 
over first few oscillations. In this range Q ai the numerator 
in Eq. (37) is still small and grows linearly with magnetic 
field. Averaging over a larger field-range similar to Ref. 
9 could not be performed because of the phase slips. 

Another interesting question is a difference in data 9 ' 11 
for the relative magnitudes G2/G1 of the second hc/2e 
and main harmonic hc/e, see Eq. (36). In the nonlinear 
transport regime this harmonic is small compared to its 
contribution in the linear transport, G2/G1 <C G2/G1, 9 
while in Ref. 11 they were comparable, G2/G1 ~ G2/G1. 
Our model also predicts the mesoscopically averaged con- 
tribution of hc/2e into linear and nonlinear conductance 
to be comparable with that of hc/e. Our approach as- 
sumes full quantum coherence of the ring, and probably 
the difference in data is due to decohcrence. 



VII. CONCLUSIONS 

In this paper, we consider mesoscopic chaotic samples 
(quantum dots or rings) and find the statistics of their 
nonlinear conductance G- This transport coefficient char- 
acterizes nonlinear DC current due to DC-bias or a recti- 
fied current due to AC bias or photon-assisted transport. 
For chaotic samples, the nonlinear effect is of quantum 
origin, which is clear from the fact that its ensemble aver- 
age over similar samples vanishes. The linear response of 
the sample in two-terminal measurements is always sym- 
metric with respect to magnetic field inversion. However, 
the Coulomb interactions lead to magnetic field asymme- 
try of the nonlinear DC response, which fluctuates due to 
the electronic interference. For the quantum dots we con- 
sider the fluctuations of (anti) symmetrized components 
Ga, G s of the nonlinear conductance. In chaotic rings the 
statistics of the phase of AB oscillations in the nonlinear 
transport regime, closely related to the ratio G a /G s , is of 
interest. 

Unlike the linear conductance measurements, in meso- 
scopic nonlinear transport experiments the way voltages 
are varied ("bias mode") turns out to be important, espe- 
cially for a weakly interacting sample. We demonstrate 
this fact qualitatively and discuss the role of Coulomb 
interactions. Quantitative self-consistent treatment of 
interactions allows us to consider magnetic-field asym- 
metry in chaotic quantum dots with many channels. Us- 
ing Eqs. (23)-(25) we show that the fluctuations of Gs 
are strongly affected by the geometry of the setup and 
discuss how the bias mode influences data of recent DC 
experiments. 

Another important issue is rectification of AC bias, 
which is quadratic in applied voltage, random, and asym- 
metric with respect to the magnetic flux inversion. The 



photovoltaic DC current can be due to rectification of 
external perturbations or quantum pumping by internal 
perturbations. Both rectification and quantum pumping 
share the aforementioned properties, and it is important 
to clearly separate them especially when the frequency of 
perturbations is high (nonadiabatic) . We consider here 
only the effects of the external perturbations and discuss 
the dependence of the fluctuations of Ga, Gs on frequency 
u>. We show that the fluctuations of both Ga and Gs, 
presented in Eqs. (27)-(30), decrease monotonically as 
uj — > 00. However, contrary to naive expectations, their 
asymptotical behavior can be very different. Since at 
high frequencies the response of the dot to the external 
bias becomes rather capacitive then resistive, the cou- 
pling to the nearby gates can be strongly enhanced. If 
the experiment is performed in a bias mode where such 
coupling contributes, the symmetrized Gl{<-°) oc 1/ui 3 
can become much larger then Ga(<-o) oc l/u 6 valid for a 
strongly interacting quantum dot. The same conclusion 
holds in the weakly interacting limit, when Q\ cx l/o> 3 / 2 
and Gl oc 1/u 3 . 

In addition, we show that recent experiments in chaotic 
Aharonov-Bohm rings might be considered similarly to 
quantum dots. The multiply connected geometry alone 
leads to AB oscillations, yet the mesoscopic distribution 
of their phase is expected to be qualitatively similar to 
that of arctan£/ /C/ s in quantum dots. Therefore, the 
bias mode should strongly affect the shape of mesoscopic 
distribution of the AB phase. The model of an AB ring, 
which we develop, consists of a dot and a long ballistic 
arm and takes into account both chaos and a ring geom- 
etry As an application of our model we consider fluctua- 
tions of the AB phase. Unlike the AB phase in the linear 
conductance, pinned to 0(mod)7r by the Onsager sym- 
metry relations, the fluctuations of the AB phase in non- 
linear transport are shown to depend on the bias mode, 
interaction strength, and temperature. 
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APPENDIX A: DIFFUSON AND COOPERON 
FOR CHAOTIC RING 

In this Appendix we determine the diffuson and 
cooperon contributions to the 5-matrix correlators of the 
random scattering matrix of a chaotic Aharonov-Bohm 
(AB) ring. This calculation is performed using Random 
Matrix Theory (RMT). 
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FIG. 7: Chaotic dot combined with long ballistic multichannel 
arm. 



First wc explain what approximations should be made 
to ensure validity of RMT. Our starting point is the as- 
sumption that the 5-matrix of the ring is uniformly dis- 
tributed over the unitary group. This means that the 
ring is essentially zero-dimensional, similarly to quantum 
dots. RMT is applicable if all energy-scales arc much 
smaller then the Thouless energy Et\i and the total flux 
through the annulus of the ring is much smaller then $o. 
Assume the ring of radius R and width W <C R to be 
diffusive with diffusion coefficient D = Ivf/2. To eval- 
uate -Exh we neglect with transversal motion of an elec- 
tron and find £rh = V r erg — {?iIvf)/2R 2 as a solution 
to Laplace equation along the circumference of the ring. 
RMT can be applied to a closed ring if the dimension- 
less conductance is large, g = i?Th/A = kplW/2R 3> 1, 
which is usually satisfied for a weak disorder even if 
W<^R. 

An open ring with ballistic contacts of the width W c 
gains a new energy parameter, the escape rate H/t^ — 
NA/2ir, where TV is the total number of ballistic chan- 
nels. The scattering matrix S is uniformly distributed 
and independent of the exact positions of the contacts 
(and therefore the length of the arms) if H/t<i -C £rh => 
ir 2 lW » 2RW C . In this case the main drop of the poten- 
tial occurs in the contacts. 

If a magnetic field is applied, the RMT is valid if the 
total flux through the annulus of the ring is much less 
then the flux quantum, <3> <C &o- Due to narrow contacts 
the time-reversal symmetry (TRS) of the 5-matrix can be 
broken at a much smaller scale, $ ~ $o \Aerg/7"d- Since 
in our rings T crg <C Td , a full crossover to the broken TRS 
can be considered. 

How well are these conditions fulfilled in the experi- 
ment? In Ref. 9 chaos was mainly due to diffusive scat- 
tering on the boundary and I w R. The width of the 
arm is 2-4 channels, while the number of channels in the 
contacts is N <~ 2, estimated from the linear conductance 
measurements, so Td/r crg ~ 5 -j- 10 » 1. In semiballis- 
tic samples of Ref. 11 (obtained by etching, and therefore 
having diffusive boundary scattering) W — W c and the 
mean free path is estimated I ~ l-^2^m^ L = 1.2/zm, the 
side length. Therefore, we have a similar estimate for the 
ratio Td/r org . Although this ratio is not parametrically 



large due to, e.g., weak disorder k?l 1, we believe that 
such AB rings still can be assumed zero-dimensional due 
to their good conducting properties together with rela- 
tively narrow contacts. 

In our calculations wc make a further simplification by 
spatially separating chaotic scattering which randomizes 
the electronic phase and the long ballistic arm attached 
to it. To find the correlators of the S- matrix elements 
we use a simplified model, see Fig. 7, which combines 
chaos and a ring geometry. A chaotic M-channel dot is 
attached to a long multi-channel ballistic arm with (M — 
N)/2 orbital channels. We assume that the size of the 
dot L and the length of the arm L a are such that L a 3> 
L 3> \J (M — N)XpL a to ensure that in the hierarchy of 
different fluxes the main flux $^ is concentrated in the 
region embraced by the arm, the flux through the dot 
$d is much smaller, but still much larger then the flux 
through the cross section of the arm. The amplitude 
of AB oscillations depends on the width of the arm oc 
(M — N). The wider the arm (relatively to the contacts) 
the closer the results should be to a uniformly chaotic 
ring. For the case when M » JV we expect it to be valid 
for the chaos uniformly distributed over the ring. Indeed, 
in this case an electron makes ~ M/N 3> 1 windings 
around the arm before exiting. 

In this appendix it is more convenient for us to work 
with an energy-dependent matrix S{e), and the final 
transformation to time-representation is rather obvious. 
The total scattering matrix S is of size N x N due to 
scattering channels in the contacts 1 and 2. Chaotic scat- 
tering in the M-channel quantum dot is characterized by 
the M x M matrix U. The scattered electron can either 
exit the sample through the N = N\ + N 2 channels (pro- 
jection operator Po = li © I2) or propagate into the arm 
with ^3 = ^4 = (M - N)/2 channels. Electrons propa- 
gate through this arm ballistically and gain phases which 
depend on the flux through the hole. In the absence of 
backscattering the electronic amplitudes at energy e are 
related to the path length L a and magnetic field phase 

o-—(; e ;)C) 

The scattering matrix of the arm is V(e) = Oi O2 P- 
Each time an electrons enters the arm cither through 
the third or fourth lead, the matrix V contributes to the 
scattering amplitude of the process. The total scattering 
matrix S is determined from the following equation: 

00 

S = Po^W'Po = PoU^-^Po, (A2) 

n=0 

where multiple n > windings are taken into account. 
Both U(e,B) and V(e 7 B) are field and energy depen- 
dent. Once we are interested only in pair correlators of 
S(e), St(e') for N, M > 1, the diffuson V and coopcron 
C of our scattering matrix are expressed via correlators 
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of the dot, T> U ,C U , and tr P(e)V* m (e'). The correlators 
of the ^/-matrix are known, see Eq. (20) for their time 
representation, and for V and C we derive 



= M -2m 



e — e 



hv F l (& d ±&; 
L 2 A 



2$n 



(A3) 



4 (A4) 



The flux penetrating the dot is denoted as $d and the 
phase <p ~ 27r$^/$o gained in the arm depends on the 
flux &h through the hole. The traces read 



tr 



x e 



(M - N) cos 

iL a [k(e)-k(s')\ 



(A5) 



Since we assumed that the area of the arm is small com- 
pared to that of the dot, the energy-dependence of Eq. 



(A5) can be neglected compared to that of T>u,Cu in 
Eq. (A4). We also assumed that since the arm is much 
longer than the size of the dot, L a 3> L, the phases <p, <fi' 
of open trajectories in the arm correspond to the flux 
§h,&'h through the hole. Therefore, the effective num- 
ber of channels Nq,v, similar to Eq. (21) for quantum 
dots is 



N c 



M - (M - N) cos 
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hv F l ( ®d±& d 
L 2 A \ 2$ 



(A6) 



The energy-dependent cooperon and diffuson in energy 
representation are given by X(e,e') = l/[Nx — 2-ni(e — 
e')/A],X = C, V. Notice that when $ = <&' the cooperon 
C is nonperiodic in the total flux <I> = & h +&d due to finite 
flux through the material of the sample, &d ^ 0. 
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